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Abstract 

The most important problem of fundamental Physics is the quantization of the gravitational 
field. A main difficulty is the lack of available experimental tests that discriminate among the 
theories proposed to quantize gravity. Recently we showed that the Standard Model(SM) itself 
contains tiny Lorentz invariance violation(LIV) terms coming from QG. All terms depend on one 
arbitrary parameter a that set the scale of QG effects. In this paper we obtain the LIV for mesons 
and nucleons and apply it to study several effects, including the GZK anomaly. 



I. INTRODUCTION 



In recent years several proposal have been advanced to select theories and predict new 
phenomena associated to the Quantum gravitational field Most of the new 

phenomenology is associated to some sort of Lorentz invariance violations(LIV's) 
Recently this approach has been subjected to severe criticism. 

In a previous letter jjj, we asserted that the main effect of QG is to deform the measure of 
integration of Feynman graphs at large four momenta by a tiny LIV. The classical lagrangian 
is unchanged. In a similar manner, we can say that QG deforms the metric of space-time, 
introducing a tiny LIV proportional to (d-4)a, d being the dimension of space time in 
Dimensional Regularization and a is the only arbitrary parameter in the model. Such 
small LIV could be due to quantum fluctuations of the metric of space-time produced by 
QG:virtual black holes as suggested inU L D-branes as in [ljj], compactification of extra- 
dimensions or spin-foam anisotropics [11 1. A precise derivation of a will have to wait for 
additional progress in the available theories of QG 291 

n 

It is possible to have modified dispersion relations without a preferred frame(DSR)[12]. 
Notice, however, that in our case the classical lagrangian is invariant under usual linear 
Lorentz transformations but not under DSR. So our LIV is more akin to radiative breaking 
of usual Lorentz symmetry than to DSR. Moreover the regulator R defined below and the 
deformed metric (5) are given in a particular inertial frame, where spatial rotational sym- 
metry is preserved. That is why, in this paper we are ascribing to the point of view of [f| 
which is widely used in the literature. The preferred frame is the one where the Cosmic 
Background Radiation is isotropic. 

Within the Standard Model, such LIV implies several remarkable effects, which are wholly 



2 



determined up to one arbitrary parameter (a). The main effects are: 

The maximal attainable velocity for particles is not the speed of light, but depends on 
the specific couplings of the particles within the Standard Model. Noticeably, this LIV 
of the dispersion relations is the only acceptable, according to the very stringent bounds 
coming from the Ultra High Energy Cosmic Rays (UHECR) spectrum [131]. Moreover, the 
specific interactions between particles in the SM, determine different maximum attainable 
velocities for each particle, a neces sary ,Zatsepin and 

Kuz'minjl^](GZK) anomaly [^l. IlI llq . Since the Auger 17] experiment is expected to pro- 
duce results in the near future, powerful tests of Lorentz invariance using the spectrum of 
UHECR will be available. 

Also birrefringence occurs for charged leptons, but not for gauge bosons. In particular, 
photons and neutrinos have different maximum attainable velocities. This could be tested 



in the next generation of neutrino detectors such as NUBE 
Vertices in the SM will pick up a finite LIV. 

This paper is organized as follows: In chapter 2, we present the LIV cutoff regulator; sec- 
tion 3 contains the effect of the regulator on One Particle Irreducible Green functions(lPI); 
section 4 defines LIV dimensional regularization ; Explicit one loop computations are con- 
tained in section 5; the LIV for mesons and baryons is found in chapter 6; Reactions thresh- 
olds are contained in section 7; Bounds on a are derived in section 8; Section 9 contains our 
conclusions. 



II. CUTOFF REGULATOR 

To see what are the implications of the asymmetry in the measure for renormalizable 
theories, we will mimic the Lorentz asymmetry of the measure by the replacement 

J d d k- > J d d kR( k2 + A f 2 ° ) 

Here R is an arbitrary function, A is a cutoff with mass dimensions, that will go to infinity at 
the end of the calculation. We normalize -R(O) = 1 to recover the original integral. R(oo) = 
to regulate the integral, a is a real parameter. Notice that we are assuming that rotational 
invariance in space is preserved. More general possibilities such as violation of rotational 
symmetry in space can be easily incorporated in our formalism. 
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This regulator has the property that for logarithmically divergent integrals, the diver- 
gent term is Lorentz invariant whereas when the cutoff goes to infinity a finite LIV part 
proportional to a remains. 



III. ONE LOOP 

Let D be the naive degree of divergence of a One Particle Irreducible (1PI) graph. The 
change in the measure induces modifications to the primitively log divergent integrals(D=0) 
In this case, the correction amounts to a finite LIV. The finite part of 1PI Green functions will 
not be affected. Therefore, Standard Model predictions are intact, except for the maximum 
arable ve,o Clty for particlesfl aad ia^ioa vertices, waica receive a finite wholly 
determined contribution from Quantum Gravity. 

Let us analyze the primitivily divergent 1PI graphs for bosons first. 

Self energy: xip) — x(0) + ^ v VtiVv + convergent, A^ v = \d^d v x{fy- We have: 



A» v = c 2 n fiu + a 

C2 is the log divergent wave function renormalization counterterm; a^ u is a finite LIV. The 
on-shell condition is: 

p 2 — m 2 — a^p^Pv = 

If spatial rotational invariance is preserved, the nonzero components of the matrix a are: 

„oo _ . a _ 

(i — cio, a — — (i\ 

So the maximum attainable velocity for this particle will be: 



c m = \i\ — — ~ 1 - (ai - a )/2 (1) 
V 1 — a 

For fermions, we have the self energy graph 

E(p) = E(0) + sT lvVvL 

s» v lv = 9 M E(0). Moreover 

s"" = sn^ + 

s is a log divergent wave function renormalization counterterm; a^ u is a finite LIV. The 
maximum attainable velocity of this particle will be given again by equation (1). 



By doing explicit computations for all particles in the SM, we get definite predictions for 
the LIV, assuming a particular regulator R. However, the dependence on R amounts to a 
multiplicative factor. So ratios of LIV's are uniquely determined. 

Vertex correction This graph has D = 0, so the regulator R will induce a tiny LIV. 

Gauge Bosons Consider the most general quadratic Lagrangian which is gauge invariant, 
but could permit LIV's |3(| 

L = c^F^F aP 

c /^o/3 j s antisymmetric in \xv and a(3 and symmetric by (a, (3) < — >(//, v) It implies that 
the most general expression for the self-energy of the gauge boson will be 

W\ V ) = (f^ Pa p^l{p) (2) 

We see that 

Pv U^{p) = 

c ^a/3 - g gj ven by a logarithmically divergent integral. We get: 

(J™*P = c 2 (r]^ a r] u/3 - V") + (3) 

C2 is a Lorentz invariant counterterm and a^ ua/3 is a LIV. 

It is clear that the same argument applies to massive gauge bosons that got their mass 
by spontaneous gauge symmetry breaking as well as to the graviton in linearized gravity. 

Explicit computations are simplified by using Dimensional Regularization as explained 
below. 



IV. LIV DIMENSIONAL REGULARIZATION 



We generalize dimensional regularization to a d dimensional space with an arbitrary 
constant metric g^. We work with a positive definite metric first and then Wick rotate. We 
will illustrate the procedure with an example. Here g = det^^) and A > 0. 



1 r d k k^ky 
^ J (2ir) d (k 2 + A) n 

y/gT(n) Jo J (27r) d m 
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1 g^T(n-l-d/2) I 



(4vr) d / 2 2 r(n) 

In the same manner, after Wick rotation, we obtain Appendix A4 of |18| |. 
These definitions preserve gauge invariance, because the integration measure is invariant 
under shifts. To get a LIV measure, we assume that 

g ^ = + (A7T) 2 ar]^r] u0 Res e=0 (5) 

where e = 2 — | and Res e=0 is the residue of the pole at e = 0. A formerly divergent integral 
will have a pole at e = 0, so when we take the physical limit, e— > 0, the answer will contain 
a LIV term. 

That is, LIV dimensional regularization consists in: 

1) Calculating the d-dimensional integrals using a general metric g^ u . 

2) Gamma matrix algebra is generalized to a general metric g^. 

3) At the end of the calculation, replace g^ v = r]^ u + (<iiT) 2 otr] tia 'rf Res e= Q and then take 
the limit e— > 0. 

To define the counterterms, we used the minimal substraction scheme(MSS); that is we 
substract the poles in e from the 1PI graphs. 

As a concrete example, let us evaluate a typical one loop integral that appears in the 
calculation of self energy graphs: 

J (2n) d [k 2 - m 2 + zO] 3 1 ; 

i <rr(2-f) i 



(47r) d / 2 2 2 ( m 2)2-| 
i rf + ^7i) 2 a5^Res €=0 T(2 - |) 1 

(47T) rf / 2 2 2 ( m 2)2-| 



2 ' ' (7) 

(8) 
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— + (ATtyad^l) + a finite LI term (9) 



4(4vr) 2V e 

LIV Dimensional Regularization reinforces our claim that these tiny LIV's originates in 
Quantum Gravity. In fact the sole change of the metric of space time is a correction of order 
e to the Minkowsky metric and this is the source of the effects studied above. Quantum 
Gravity is the strongest candidate to produce such effects because the gravitational field is 
precisely the metric of space-time and tiny LIV modifications to the flat Minkowsky metric 
may be produced by quantum fluctuations. 
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V. EXPLICIT ONE LOOP COMPUTATIONS 



A. Example:^ 3 in six space-time dimensions 

To illustrate the method using the cutoff regulator, we consider g<ft 3 in six space-time 
dimensions. 

The one-loop contribution to the self energy of the particle is: 

{-gif f d 6 k 1 1 

mq) - — ^— J ( 27r )6 k 2 _ m 2 + i0 _ q y _ m 2 + i0 y 10 > 

The term containing the LIV is: 



L(i X ) = -2q^g 2 B^ (11) 

r> f d 6 k k^k v 

^ J (2tt) 6 (k 2 - m 2 + iO) 4 1 ' 

To evaluate B^, introduce the regulator of the integration measure, 

-A 2 

R= k 2_ A 2 + ak 2 + lQ ( 13 ) 

define k = Ap and take the limit A— > oo. In this way we verify that the LIV is mass 
independent. Since a « 1, we keep only the first order in a. We end up with: 

r R „ f d6 P PoP&v /i.X 

^ a J (2tt)6 (p2 - 1 + * ) 2 (p 2 + *0) 4 1 } 



Therefore: 



g 2 aq 2 i 

L[%X) = (15) 



B. LIV in the Standard Model 



We follow 20|, |2l( and use LIV Dimensional Regularization. 
Photons 

In the SM the photon self-energy can be written: 



-T\e 2 



so that the LIV photon self-energy in the SM is: 



23 

LW u {q) = -—e 2 aq a q p 



(r^« + rTSSSg - if08S8S - vT& v M (17) 
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It follows that the maximal attainable velocity is 

23 

c 7 = 1 - ye 2 a (18) 

We have included coupling to quarks and charged leptons as well as 3 generations and 
color. 



C. Fermions 

Let us consider QED, as an example. The electron self-energy to one loop is given by: 

2 r d d k 1 (iU + m) —i . . 

- ,E 2 «z) = He) J J^- d ^r k2 _ m2 + iQ l, {k _ q)2 _ fx2 + i0 (19) 

To obtain the LIV, we have to evaluate (we have introduced a parameter A and put it to 
zero afterwards): 

- iLE 2 (q) = 2i(-ie) 2 / - — -r — , ' - NQ 20 

iyi> v ; J 2tt V* /o k 2 - A + iO 3 V ; 



(2ir) d y/g (k 2 - A + iO) 3 

= -^ d -^S§fT9W^W (21) 



2 (rf-2) 1 r(2-d/2) 



but 



^ 2 (4*)^ ' (22) 



f r ( 2 ~ rf / 2 ) _ a ,p _ M 2 aRes £=0 T(2-d/2) _ 

4- ^5.7° (24) 



SO, 



2 

-^LS 2 (g) = ^go7° (25) 
Similarly, in the SM, the fermion self-energy is given by: 

(47r 2 )S(g) = - -4 (\c v + c A \ 2 P L + (\c v -c A \ 2 P R ) + finite (26) 

^ graphs 

where the fermion-gauge boson vertex is: 

*7 fc (cv-CA7 5 ) (27) 
and P L (P R ) are the L(R) helicity projectors. 
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Therefore 

LTJn) = 



LE(g) = % ol ° ]T (\ c v + ca\ 2 Pl + i\c v - c A \ 2 P R ) (28) 



graphs 

We apply this last result to neutrinos and charged leptons below. 
Neutrinos: The maximal attainable velocity is 

c u = l-(3 + tan 2 9 w )^— (29) 

8 

In this scenario, we predict that neutrinos [ljjj emitted simultaneously with photons in 
gamma ray bursts will not arrive simultaneously to Earth . The time delay during a flight 
from a source situated at a distance D will be of the order of (5 x 10 _23 )-D/c ~ 5 x 10~ 6 
s, assuming D = 10 10 light-years. No dependence of the time delay on the energy of high 
energy photons or neutrinos should be observed(contrast with [1]). Photons will arrive 
earlier since a < (See below). These predictions could be tested in the next generation of 
neutrino detectors 

Using i?£-gauges we have checked that the LIV is gauge invariant. The gauge parameter 
affects the Lorentz invariant part only. 

Electron self-energy in the Weinberg-Salam model. Birrefringence: 

Define: ei = ^-e, e# = ^y-e, where e is the electron field. We get 

c L = l- {-^{sin 2 9 w - 1/2) 2 + e 2 + g 2 /2)^ (30) 
cos z w 2 

CR = 1 ' {e + ^e7 ) 2 (31) 

The difference in maximal speed for the left and right helicities is ~ (5 x 10~ 24 ). 

VI. MESONS AND BARYONS 

In order to apply our results to the computation of the UHECR spectrum and other 
phenomena, we must calculate the maximal attainable velocity of hadrons. As we mentioned 
before, the problem is hadronization. One way to get an estimation of the effect is using 
effective lagrangians. 

We US e the rcSU , tS of H Q fer the wave functl on —Nation of pl on S and oucleons 
in the chiral lagrangian and Heavy Baryon Chiral Perturbation Theory. They get: 
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* 1 = 1 -W$F*l + »*' (32) 

Here, m n is the renormalized pion mass, F is the renormalized decay constant of pions 
and gA is the axial vector coupling constant, in the chiral limit. 

Using the LIV metric, we can read off the maximal attainable velocities for pions and 
nucleons: 

2m!a 



C,r = 1 + 



3F 2 



c„ = l + ^ (34) 



VII. REACTION THRESHOLDS 

Knowing the LIV for nucleons, pions, photons and electrons, we proceed to study the 
reactions involved in the GZK cutoff. We follow the discussion in jflQ]. 

A. Photo-Pion Production 7 + p —* p + ir 

Let us begin with the photo-pion production 7 + p — > p + tt. Considering the corrections 
provided in the dispersion relation (J34j) for pions and nucleons, we note that, for the photo- 
pion production to proceed, the following condition must be satisfied 

2 5cE% + 4E^ > "*( 2 ™p + ™») . (35) 
where E n is the produced pion energy and 5c = c p — c w . 

B. Pair Creation j + p — > p + e + + e~ 

Pair creation, 7 + p — > p + e + + e~, is greatly abundant in the sector previous to the 
GZK limit. When the dispersion relations for fermions are considered for both protons and 
electrons, it is possible to find 

5c — -E 2 + Eoj > m e (m v + m e ), (36) 
m p 

where E is the incident proton energy and 5c = c p — c e . 
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VIII. BOUNDS 



In order to study the threshold conditions (1351) and (IHoT) . in the context of the GZK 

nn 

anomaly, we must establish some criteria. Firstly, as it is seen in [13, 24], the conventionally 
obtained theoretical spectrum provides a very good description of the phenomena up to an 
energy ~ 4 x 10 19 eV. The main reaction taking place in this well described region is the pair 
creation 7 + p — » p + e + + e~ and, therefore, no modifications are present for this reaction up 
to ~ 4 x 10 19 eV. As a consequence, and since threshold conditions offer a measure of how 
modified kinematics is, we will require that the threshold condition (|3*Bj) for pair creation 
not be substantially altered by the new corrective terms. 

Secondly, we have the GZK anomaly itself, which we want to explain. Since for energies 
greater than ~ 8 x 10 19 eV the conventional theoretical spectrum does not fit the exper- 
imental data well, we shall require that QG corrections be able to offer a violation of the 
GZK-cutoff. The dominant reaction in the violated E > 8 x 10 19 region is the photo-pion 
production and, therefore, we further require that the new corrective terms present in the 
kinematical calculations be able to shift the threshold significantly to preclude the reaction. 

We begin our analysis with the threshold condition for pair production. In this case we 
have: 

5c — -E 2 + Ecu > m e (m v + m e ), (37) 
m p 

with 5c = c p — c e . As is clear from the above condition, the minimum soft-photon energy 
co> min for the pair production to occur, is 

^min = -ttO^p + rn e ) - 5c —E. (38) 
h/ Trip 

It follows therefore that the condition for a significant increase or decrease in the threshold 
energy for pair production becomes \5c\ > mp{m p + m e )/E 2 . In this way, if we do not want 
kinematics to be modified up to a reference energy E re f = 3 x 10 19 , we must impose the 
following constraint 



|c p -c e |< K + r e)mp =9.8xlO- 22 . (39) 



Similar treatments can be found for the analysis of other astrophysical signals like the Mkn 
501 7-rays j^j, when the absence of anomalies is considered. 
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Let us now consider the threshold condition for the photo-pion production. We have 

2 6c El + AE n u > — ^ p ; . 40 

m p + m n 

It is possible to find that for the above condition to be violated for all energies E n of the 
emerging pion, and therefore no reaction to take place, the following inequality must hold 

2a; 2 (m„ + m n ) 0A r . l9 . 

C - - °P > 2fO i ^ = 3 ' 3 X 10 ^ • 41 

m^2m p + m ?r ) 

where wo = KT = 2.35 x 10~ 4 eV is the thermal CMBR energy. 

Combining the two reactions and the standard values, = 139Mev, = 1.26, F = 
92AMev, we get an upper and lower bound on a 



2.2 x 1(T 21 > -a > 1.3 x 1(T 24 (42) 

First of all, we notice that a < 0, in order to suppress the photopion production, thus 
removing the GZK cutoff. This implies that photons are the fastest particles and they arrive 
before neutrinos coming from the same source of GRB. Moreover, photons become unstable. 

n 

They decay in a electron positron pair above an energy £7oH- See below. 

Since c p h oton > c prot0 n, the strong bound of (2?! is avoided: Proton is stable under Cerenkov 
radiation in vacuum. 

If no GZK anomaly is confirmed in future experimental observations, then we should 
state a stronger bound for the difference c n — c p . Using the same assumptions to set the 
restriction (|39J) when the primordial proton reference energy is E Te f = 2 x 10 20 eV, it is 
possible to find 

\cn - c p \ < 2.3 x 1(T 23 . (43) 

In terms of a, this last bound may be read as 

\a\ < 9.1 x 1(T 24 , (44) 

which is a stronger bound over a than (|39p. offered by pair creation. 
Photon unstability 

It has been pointed out in (f| Q that if c p hoton > c e i ectr on then the process 7 — ► e + + e~ 



is allowed above an energy Eq\ 



/ 2 

E = m e J— (45) 
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where 5c = c 7 — c e . 



In our case, we have: 

6 cos 2 9 



23 r7 2 

5c L = -«(-e 2 - (-^— (sinX - 1/2) 2 + e 2 + //2)/2) (46) 



** - -« ( f 62 - + "^t^ (47) 



Therefore, with 



EL = 2.3 x 10 8 Get; 

ER = 1.9 x 10 8 Gew (48) 

So, we should not detect photons with energies above 2.3 x 10 8 Get> 
Neutral pion Stability 

Following we study the main decay process of neutral pion 7r — > 7 + 7 .This becomes 
possible if c 7 > and above an energy 

E n = , (49) 
2(c 7 c^) 



Using the bound c 7 — c w < 10 22 obtained in 



28j | . we get 



|a| < 5.4 x 10" 23 ( 50) 

In our numerical estimates we have chosen a = —5 x 10 -23 . 

We get E n = 10 19 e\A Therefore we expect that neutral pions above this energy are 
stable, so they could be a primary component of UHECR. Photons will be unstable above 
this energy by the same mechanism. Notice however that photons are unstable at a lower 
energy due to electron-positron pair creation 



IX. CONCLUSIONS 



In this paper we have computed the LIV induced by Quantum Gravity on Baryons and 
Mesons, using the Chiral Lagrangian approach. This permitted to fix that a < 0, in order 
to explain the GZK anomaly. Studying several available processes, we found bounds on a: 

From pair creation and absence of photopion creation: 2.2 x 10~ 21 > —a > 1.3 x 10 -24 . 

From pion stability and the most stringent experimental bound found in |28j: |a| < 
5.4 x 10~ 23 . 
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Then, several predictions are obtained:Photons are unstable above an energy 2.3 x W 8 Gev . 

Neutral pions are stable above an energy E n = 10 19 eV; so they could be a primary 
component of UHECR, thus evading the GZK cutoff. 

Moreover, in time of flight experiments, photons will arrive before neutrinos, assuming 
that they were emitted simultaneously at the source. No energy dependence of the time delay 
should be observed. The time delay during a flight from a source situated at a distance D 
will be of the order of (5 x lO~ 23 )D/c ~ 5 x 10 -6 s, assuming D = 10 10 light-years. 
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